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Abstract 

We calculate the decay rates of leading Regge trajectory states 
for very high level number in open bosonic string theories, ignoring 
tachyon final states. The optical theorem simplifies the analysis while 
enabling identification of the different mass level decay channels. Our 
main result is that (in four dimensions) the greatest single channel 
is the emission of a single photon and a state of the next mass level 
down. A simple asymptotic formula for arbitrarily high level number is 
given for this process. We also calculate the total decay rate exactly 
up to N=100. It shows little variation over this range but appears 
to decrease for larger N. We check our formalism in examples and 
calculate the decay rate of the first excited level for open superstring 
theories. The calculation may also have implications for high spin 
meson resonances. 

1 On leave from: Theoretical Physics Group, Blackett Laboratory, Imperial College, 
London S.W. 7 2BZ; After September 1st: Joseph Henry Laboratory, Princeton University, 
Princeton, NJ 08544. 

2 After September 1st: Joseph Henry Laboratory, Princeton University, Princeton, NJ 
08544. 

3 On leave from: Theoretical Physics Group, Blackett Laboratory, Imperial College, 
London S.W. 7 2BZ. 

4 After October 1st: CERN Theory Division, 1211 Geneva 23, Switzerland. 


O 

CM 


>• H 
.x: a 


a o 





LJ 

•H 

<5 M 

O CD 

M — 

iJU LQ 
tr* O 


^ H ^ 

s/J ' U 

co J 

i O -t* 

03 jq M 
f OJ o 

I o 

j u u 
I *1 W 

V) H CD 

u u 


Operated by Universities Research Association Inc. under contract with the United States Department of Energy 


Introduction 

In all the years of string theory [1], remarkably little attention has been 
paid to the interactions of the massive string modes. Recent work [2,3,4] 
on the statistical mechanics of strings at very high densities has emphasised 
the fact that as the Hagedorn temperature is approached, a phase transition 
occurs with most of the energy going into very long strings. It is therefore of 
some importance to understand the interactions of very excited strings and 
in particular how fast they decay. 

Quite independently of this it is interesting to try and understand whether 
and how a classical limit emerges in string theory - do highly excited strings 
behave in a ‘classical 5 way? 

In this paper we attempt a modest beginning to the effort at a better 
understanding of the massive string states by calculating the decay rate of 
a particular class of occupation number states of arbitrarily high energy - 
the states on the leading Regge trajectory of open string theories. A similar 
problem was considered in 1971 by Green and Veneziano [5] who argued 
that the high mass intermediate states in dual amplitudes had narrow decay 
widths and our result appears consistent with theirs. However their argument 
was rather incomplete. 

Our main tool will be the optical theorem and a method for extracting 
the imaginary part of the string self-energy diagram which we explain in the 
first section. One advantage of our method is that it enables one to identify 
which piece of the imaginary part is responsible for the decay into a given 
mass level for each of the final states. Thus while we are saved considerable 
labour in calculating polarization sums and phase space integrals we can 
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recover quite a lot of detailed information on how massive strings decay. 

This is actually essential to how we treat the tachyon. Of course tachyons, 
living on a continuous Lorentz hyperboloid 1 , p 2 = m 2 , do not have a definite 
sign for their energy. Any decay process involving tachyons is infinite - one 
can obtain in the decay an arbitrarily negative energy tachyon plus an arbi- 
trarily positive energy excited state in the ‘decay’ of a photon for example. 
This is quite a graphic way of seeing the ‘vacuum’ instability in theories with 
tachyons. 

However we wish nevertheless to calculate physically meaningful quanti- 
ties for bosonic strings since these are certainly the most tractable and allow 
us to proceed furthest. Since these theories have tachyons we will simply 
calculate the total decay rate of highly excited strings into everything ex- 
cept tachyons. A major advantage of our method is that it enables us to 
extract the ‘tachyonic’ part of the decay rate rather simply and calculate the 
remaining, finite part. This is perfectly well defined and in fact it is quite 
likely that the main features of our calculation will remain true for theories 
without tachyons - heterotic strings and superstrings for example. 

We shall calculate in arbitrary d. The open bosonic string is of course 
inconsistent at the one loop level outside of d = 26, and in fact the need for 
d — 26 was discovered in calculations by Lovelace[6] of one of the diagrams 
we consider. However we shall really only be using the loop diagrams as a 
convenient way of summing tree diagrams (we shall only use it s imaginary 
part and not the real part which is infinite), which do not require d = 26. 
Furthermore, because of the particular external states we consider the loop 

1 Our metric convertion throughout this paper is (-,+)+,+•••)• 


3 


ORIGINAL PAGE IS 
OF POOR QUALITY 


amplitude will not have the singularities of the type considered by Lovelace, 
corresponding to propagating closed strings. Our calculations would also 
apply to the case of open strings in d = 26 with a subset of the 26 dimensions 
taken to be toroidally compactified, and with the compactification radius so 
small that degrees of freedom depending on the extra dimensions are never 
excited. 

We shall discuss in detail the diagrammatics for the simplest case, the 
orientable U(l) open string. The extension of our result to other groups 
should be straightforward and not qualitatively different. 

Our main result is that the largest single process is the decay into a 
photon and a state at the next level down - for this process the decay rate 
is inversely proportional to the mass of the initial string state. The total 
decay rate is much harder to say anything analytic about. We calculate it 
up to N=100 and it appears to be slowly decreasing for large N. This means 
that the open string coupling constant cannot be thought of as a splitting 
probability per unit length - if this were the case long strings would fragment 
at a rate increasing as \^N. 

There may well be implications of our result for hadronic physics. It 
is well known that resonances of very high spin are observed rather more 
commonly than would be expected naively - their widths being generally no 
larger than resonances of lower mass. 

One would expect that the bosonic string description of mesons would be 
good for the leading trajectory resonances since these are the most spatially 
extended states. Our procedure of throwing away the tachyon is certainly a 
fairly violent procedure at the low end of the mass spectrum so our calcula- 
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tions are probably a poor description of the interactions of low mass states. 
However the decay rates we find are not dominated by the emission of the 
lowest mass states as can be seen for example in Figure 7. It seems reason- 
able therefore that our calculations should provide a reasonable decription 
of the decays of high spin mesons. It is gratifying that our results lead to an 
approximately constant decay rate as a function of level number, in qualita- 
tive agreement with experiment. It would be interesting to try to make the 
connection more precise - this could provide new information about hadronic 
strings. 

As a final motivation for this work, it may eventually provide a consistent 
method for treating ‘back-reaction’ problems in radiation from cosmic strings 
- gravitational, electromagnetic or goldstone radiation. In the quantised 
string these processes are automatically finite, and the probability amplitudes 
for various final states should also tell us the nature of the configuration that 
the string ends up in. We leave this however for future work. 

In a subsequent paper three of us [7] will extend the treatment presented 
below to closed strings where interesting issues connected with modular in- 
variance arise. 

The outline of this paper is as follows. In Section 1 we give a simple exam- 
ple of the optical theorem in d-dimensional scalar field theory and a method 
of extracting imaginary parts which we later use for strings. In Section 2 we 
calculate the one loop planar self-energy diagram for arbitrary highly excited 
external strings on the leading Regge trajectory, for which the vertex oper- 
ators are particularly simple. In Section 3 we calculate its imaginary part, 
and in section 4 we discuss how the other topologically distinct diagram at 
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this order should be added. The total decay rate as a function of mass of 
the initial state is calculated exactly up to level 100, illustrated in Figures 6 
to 8. In Section 5 we find an asymptotic formula for the single greatest term 
in the decay, due to the string emitting a photon and lowering its mass level 
by one. In Section 6 we calculate the decay of the first massive state in open 
superstring theories by the same method. 

While this work was being completed we recieved a preprint by J. Polchin- 
ski [12] who discusses the interactions of macroscopic closed strings on a 
torus. He mentions that by a similar method one could calculate the decay 
of a closed string into open strings. His result seems quite different to ours 
- he claims it is interpretable as a splitting probability per unit length. It 
would be interesting to understand the reason for this apparent discrepancy. 

1 The Optical Theorem: A Simple Example 

In this section we present a simple example of the optical theorem for inter- 
acting scalar fields. This is of course very well known and in this example 
rather unnecessary - integrating over phase space is simpler. However the 
method we use will generalise directly to the case of strings, both open (in 
this paper) and closed (in a sequel [7]). It may also seem rather pedestrian 
for open strings, where a quicker approach would be to use a generalisation 
of the Cutkosky rules (see for example [11]). However both in the interests of 
giving a rigorous self contained account and because our approach generalises 
to closed strings, we prefer the method we shall explain below. 

To recap briefly, the optical theorem starts from the unitarity of the S 
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S = 1 + iT sts = l-»i(Tt-T) = TtT (1) 

Taking matrix elements and defining < / | T j i >= T/^(p/ — p;) we obtain 2 

2lmTa = ^2 | Tif |V(p/-pO (2) 

/ 

where of course initial and final states must be normalised, producing the 
usual 2 E factors and the volume and time factors which turn the right hand 
side into a decay rate. 

Now consider a theory with two scalar fields $ and <f> with masses M 
and m respectively, and a interaction. Then equation (2) relates the 

imaginary part of the $ self-energy to the total decay rate of $ particles into 
<j> particles . To lowest order in the coupling, this is shown diagrammatically 
in Figure 1. 

In its rest frame, the $ particle has momentum k M = (ikf, 0) and we have 
for the one loop $ self-energy 

A 2 r tPp 




(p 2 + tti 2 — ie)((& — p) 2 4* m 2 — ie) 


( 3 ) 


where ^ is from normalising the $ particle initial state and | from the 
symmetry factor. 

Using Schwinger’s representation this is 
A 2 


- — f [°° dct /“dy3e" ia(pJ+mJ_ “ )_ ‘ /3((fc " p))a+m3_ ‘‘ ) 
4 M J Jo Jo 


( 4 ) 


which is of course divergent for d > 4. However from (2) we know that its 
imaginary part, being a phase space integral over tree diagrams, must be 
2 our convention throughout is tfp = and tf(p) = 2 *6(p) 
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finite. Performing a Wick rotation, p° = ip° E we regularise (4) with a Eu- 
clidean convergence factor e -Ap ^. We shall extract the imaginary part of the 
regularised expression and then set A to zero. The Schwinger representation 
arises naturally in string theory as we shall see, and this method of regular- 
isation is close to what happens in the closed string theory, in which there 
are no ultraviolet divergences. After performing the Gaussian momentum 
integral we change variables to x = a — iA + (3,y = to obtain 

T - 

A tt 


My) = 

where we drop an overall real constant e A3 which tends to 1 as A tends to 
zero. Taking the y integral along the real axis, the x integral converges for 
all A(y ) on the contour (1) shown in Figure 2. 

Now the function A(y ) is a parabola (Figure 3); for M 2 < 4m 2 it is always 
positive there is no imaginary part and consequently no decay - simply a 
reflection of energy conservation. However if M 2 > 4m 2 then A(y) is negative 
over a range of y which results in an imaginary part for the self energy and 
the $ particle decays. 

To see this is so, for A(y) positive we can rotate the x contour downward 
to (2), so x = —iz, A < z < oo. Then Tu is now clearly real ( although of 
course divergent as A tends to zero) and there is no decay. 

However if M 2 > 4m 2 then A(y) is negative between y± = |(l±\/l — 
and the contour may only be rotated upward , running around the origin on 
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m 2 - k 2 y(l - y) 


( 5 ) 
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(3). Integrating by parts, we can reduce I(y) in (4) to a real series in negative 
powers of A, and a remaining integral. 

If d is even, setting z =| A | x we find 


I = 


= I A I* 


A l^ 2 (0 
(( 


a-f 


i: 


dzz 1 ~*e iz 


2 - - 
r A 2 


|A|A 


+ 


V 2 


(2-f) (2 - f )(3 - f ) 


+ ••■)« 


+ 


1 


r(f-i)v. z 


(J d in) 


( 6 ) 


where r =| A | A. As before, the series (in inverse powers of A) is reed. 
However for finite A the integral along (3) may be distorted to run up along 
the imaginary axis, along an infinitesimal semicircle around the origin and 
along the imaginary axis up to ioo. The integral along the axis gives a 
real principal value (logarithmically divergent as A tends to zero), while the 
integral around the semicircle gives tV. Thus for A negative 

A I*" 2 . 


Im(I) = 


r(f-i) 


ITT 


(7) 


( 8 ) 


Similarly if d is odd we obtain a real series ending in 

. I A |j -2 r- f°° dx _ x 

where we have distorted the integral to run along the imaginary axis. The 
integral is now finite as A tends to zero and the result is again given by (7). 
These results are also easily obtained in dimensional regularisation. 

The remaining y integral in (5) may be performed as explained in Ap- 
pendix A to obtain 


A 4m 1 

(le^rt^r M’’ 


d - 3 
2 


(9) 
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It is simple enough to check that this is identical with the answer obtained 
from the right hand side of (2) 

A 2 1 


r = 


2 2 M 


> 1 (io > 


The | comes from the identical <f> final particles, with momentum pi and p 2 * 
One slightly surprising consequence of (9) is that decay rates decrease 
rapidly in high dimensions, if the masses and the dimensionless coupling 
A mV are kept fixed. This is a result of the fact that the surface area of a 
unit sphere decreases rapidly in high dimensions. For just this trivial reason 
massive string modes are rather long lived in 26 dimensions! 


One Loop String Amplitudes and Decay 
Rates 


We now turn to the problem of evaluating one loop string amplitudes with 
two identical excited external particles. We shall focus on open bosonic 
string theories with external particles of definite occupation number, and 
futhermore on the states with highest angular momentum for given mass, the 
“leading Regge trajectory” states. Classically they correspond to a straight 
rotating rod where of course the ends move at the speed of light, although 
they are of course occupation number rather than position eigenstates. 

We construct the one loop amplitude by inserting two vertex operators 
to produce the external states and then tracing over all states circulating 
in the loop. This is a well known procedure for external tachyons or vector 
particles. In this section we present the calculation for arbitrarily excited 
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external states and although for the most part this is mostly familiar, there 
are some new aspects and we shall highlight these. 

We work in d-dimensional spacetime - we shall not see the requirement 
d=26 in this process as we discuss later. The states we calculate the decay 
of are, in standard notation[l] at level N 

| >N= ?W..a a - I 0 > ( n ) 

with momentum k obeying the mass shell condition k 2 — 2 tcT(N 1) and 
the (symmetric) polarisation tensor 77 has N indices. T is the string tension. 
In what follows we will set TV = = 1 and restore T by dimensional 

analysis later. Requiring (11) to be physical (i.e. L n \ ij> >jv= 0,n > 0 ) 
results in the transversality and tracelessness conditions 

= 0 ( 12 ) 

The vertex operator for these states is simply 

14 = W „ T : P p P v P p ....P T e ikX : (13) 

where dots denote normal ordering and the notation is standard[l]. 

The conditions (12) are sufficient, with the mass-shell condition, to guar- 
antee (13) has the correct conformal dimension. There are no normal ordering 
singularities between any of the different terms in (13), because of (12), so 
the naive conformal dimension is the correct one. It is also easily checked 
that applying the vertex (13) to the vacuum at t = ioo in the usual way [1] 
produces the state (11). Again following standard procedure we construct 
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the one loop amplitude by sewing together the ends of the tree diagram (Fig- 
ure 4) with loop momentum p. Including ghosts in the propagators [ 9 , 1 ], we 
obtain 


iT = g 2 J dx-i J dx 2 J cPp Tr(x[°V k x 2 °V. k )x 1 2 “x 2 2_,e (14) 

The evaluation of (14) is well known for external tachyons or photons[lO,l] 
and we shall not repeat it here. Generalising this method to the case at hand 
we write 


Y _ ® \ . ik.X+l.P I 

V ■' £=0 


(15) 


This allows considerable simplification of the trace calculation. Indeed we 


now arrive at 


iT 

G 

H 


2 ^ 

= ^4 Jo dxi lo dXi J x i 2lt x 2 2 ~ U x i 12 X 2 P_fc),/2 f 2 ~ d (w) G 


A -UT...K. 

V V 


= eS”=i -'*•(*»£) 


(—....—) (-2- —) e (p e t.(p-k) H 

lf ={=0 


(16) 


where rj denotes the complex conjugate of 77 and 


/(«) 


XjX 2 


na-o 

71= 1 

+ n <{ + - *J) 

+»’«>? + *J) + (('■ ‘ + (17) 


We note that due to the conditions of transversality and tracelessness of the 
polarisation tensor, (16) can be rewritten as 


G = T}^- A 7 fT...« 


(— —)(-£- -2-)e {p e fp O I - 
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Q = e EZ L l n(*r+xJK.«7(l-«'“) 


Let us define B = Performing the differentiations and paying 

careful attention to the combinatorial factors we arrive at 


G = f (wr 


N indices 


r terms 


err % T ... b 


N-r 


(riy(N-r)l 

with all “spare” indices being contracted. Hence (16) can be written 


iT = J dxi J dx 2 J x t 


2 -it 2— «'« JP 1 !" 1 ^.(P-*) 1 / 2 f2~d 


x\ x 


r~ d {w ) u 


u = £ 

r= 0 
oo 

^ = n 


* (IV!) 2 


- Tf v - f T -(py...)(pV...) 

(r!) 2 (lV - r)! 7 7 ^ P ^ } 

(1 — Xim m-1 )(l — x 2 w m ~ 1 ) 


m=l (1 — w m ) 2 

Now we change variables to the ’Schwinger’ variables a and /? 


(18) 


Xi = € 


x 2 = e 1/3 


Next we Wick rotate and regularise with a Euclidean convergence factor 
e -Ap « as we did in Section 1. Performing the Gaussian integrals and defining 
x = a — t'A + /?, y = again as in Section 1 we obtain 


T 

a 

P 

C 


_i jy V [ l dy [°° dx (i) 2+r " d/J *-'" d/3+1 e-“*P 

2P ext 4(27r) rf / 2 Jo 47 J—iA V 7 

[ 3/( 1 -y)(l -IV) + 1 -ze] 

(W!) 


/(tn) 2 d C(xi,x 2 ) 2(7V 1} P(x 1 ,x 2 ) Af r 


r!(JV-r)! 


(ft(n+l) tt(n-l) 


fw~ 3 ( E 1 ** • ) 


(19) 
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Here we have set T] 2 = 2E 1 ^ in order to normalise the external states ( 11 ) 
and have used Jacobi’s identity [ 8 ] to reexpress the infinite product C in 
terms of the partition function and an infinite series ( a Jacobi 6 4 function). 

This is our final expression for the one loop amplitude with identical ex- 
cited external particles. We shall now determine how to extract the imaginary 
part of it just as in the field theory example. 

3 The Imaginary Part of the String One Loop 
amplitude 


Let us concentrate on the expression for the amplitude (19). In particular 
let us consider the factor P. This is simply a polynomial in Xi and the 
term having a coefficient that also depends on the summation variable 

r i.e., 

~ P, * = 0,1,2,... (20) 


p= £ *; 

p,q-0 


P T <1 
X 2 


where the C* are numbers. 

PQ 


In terms of x and y we have, recalling from 19 that a = xy + xA,/3 = 

*(i - y ), 


Xi = e~ lxv 


x 2 = 


where we drop an overzdl real constant e A in which tends to unity as A 
tends to zero. This allows us to include the factor of P in (19) as 

\ 2 +r-d /2 x -T-d/ 2+1 K 


T - 


2E, 

J( = e -iA(y)x-tx 




( 21 ) 


14 



ORIGINAL PAGE IS 
Of POOR QUALITY 


following the field theory example we set 

A(y) = —(N - l)y(l - y) - 1+ W + ?(! - v) ( 22 ) 

The positive and negative regions of this function will determine whether or 
not our string can decay. However the -1 in (22) is a reflection of the tachyon. 
if this is not removed then the decay rate will be infinite. To see this, return 
to (14) where the integrand contains 

Tr(xf° V k 4° VLfc) 

and the trace is over states propagating around the loop. Now xf 0 is just 
x f /2JrN where N is the level number operator. Inserting a complete set of 
occupation number states one sees that the power of xj occurring in the 
trace corresponds to the level number of the particle circulating on that side 
of the loop, and similarly for X 2 . In particular terms independent of X\ or X 2 
correspond to tachyonic states propagating on either side of the loop. Thus 
if we wish to extract those parts of P which do not correspond to tachyonic 
decays we simply discard the terms in P which do not have at least one power 
of X\ and x 2 . We emphasise that this results in a perfectly well defined (and 
totally finite) calculation, and corresponds to summing up all tree diagrams 
not involving tachyons and integrating over phase space. We perform an 
explicit check of this for the simplest case, the decay of the first massive 
state into photons, in Appendix B. 

The only terms in P we keep are those with at least one power of x\ and 
x 2 . This results in the -1 in (22) being cancelled. One can now perform the 


15 


x integral and take the imaginary part just as in section 1 to get 


Im(T) = 


1 g 2 7T 

2 E ext 4(2t r)*/* 


N 




r iy I iW 

4- r(f+r-l) 


( 23 ) 


For what values of p and q does A(y ) go negative?. We find, for jV > 1 
that this happens for 


(P~ g) 2 _ 2 (p + g - 2) 

(A - l) 2 (AT - 1) 


(24) 


This is the string generalisation of the constraint k 2 > 4m 2 we had before. 
The values of p and q correspond to different mass levels in the decay products 
and (24) just tells us if certain decays are energetically allowed. This is not 
the full story - for a given p, q that satisfies (24) we calculate the range of y 
for which A is negative (see Figure 3) 


y+ -y- = 


(p-q) 

(N-iy 


2(p + q — 2 ) 

(TV - 1) 


+ 1] = 2a 


(25) 


Returning to (23) we have to evaluate the y integral, just as in the field 
theory example. From Appendix A we find 





where a is defined in (25). 
Thus (23) becomes 


Im(T) = 




g i t 


2E ext 4(27r) d / 2 £ t 


jrj2c; - 1 '> hr ' 1 


PI 


r(i + r-l) 


(26) 


where the p, q sums only run over values for which the energy constraint (24) 
holds. 
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4 Twisted Diagrams 

Equation (26) is not the final answer because it comes from only one diagram, 
the planar diagram shown in Figure 4. In general for open strings there are 
two other classes of diagrams, the orientable and nonorientable nonplanar di- 
agrams. These are illustrated in Figure 5. To construct these diagrams one 
has to use the twist operator. Let us assume for simplicity that we are deal- 
ing with (orientable) U(l) strings - which amounts to imagining the strings 
have two oppositely charged particles (’quarks’) on either end . Charge con- 
jugation symmetry (C) reverses the string - it is the same as the operation 
(—1)^ where N is the level number operator. For example tachyons may be 
thought of as neutral point particles, even under C. Photons are odd under 
C and so on[l]. 

To see what this means in tree amplitudes, consider the general 3 point 
vertex. The general prescription in string theory is to add all cyclically 
inequivalent diagrams contributing to a given process. In this case there are 
two diagrams, < 1 | V(2) | 3 > and < 3 | V(2) | 1 > in bra and ket notation. 
These correspond to the string | 3 > breaking in 2 ways - emitting state 2 
from one end and state 1 from the other and the reverse. The second diagram 
is however related to the first: 

< 3 | V{2) | 1 >= (-l)*+*«+*» < 1 | V(2 ) | 3 > (27) 

by a generalisation of the cyclic symmetry proof. Here Ni is the level number 
of the ith state. 

This means that the total amplitude of a given process is zero unless 
(— 1)^ J = ( — i)Wj+Wi # an expression of C conservation. For example the 
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amplitude for a tachyon (and more generally any string state) to scatter off a 
photon is zero, indicating that the string states carry no net charge, consistent 
with the picture of their coupling to photons as “dipoles”, carrying oppositely 
charged particles at either end. 

The second diagram may be written in terms of the first diagram by using 
the twist operator S7 = ( — l) 7 ^ 

< 3 | V{2) | 1 >= (-1)*’ < 1 | ftV(2)ft | 3 > (28) 

How does this relate to loop diagrams? For oriented open strings, the nonori- 
entable diagram (Figure 5a) does not contribute - cutting the loop reveals a 
string with two like charges at the ends. So there are just the planar and 
nonplanar orientable (Figure 5b) diagrams to consider. Now it is seen that 
squaring the two tree diagrams and summing over intermediate states will 
result in two direct terms, equal to the planar diagram, and two interference 
terms, equal to the nonplanar orientable diagram. This shows that for our 
loop diagrams to give a unitary result, we must add in the loop diagram 
with two twists with the same weight as the planar diagram but with a level- 
dependent phase factor ( — 1) N as well. Thus for the nonplanar orientable 
diagram the trace in (14) is replaced with 

Tr(x^nv k x^nv_ k )(-lf (29) 

where N is the level number of the external particle. 

The result of all this is very simple. The effect of the twist operator in 
the trace part of the loop amplitude is to change Xi and x 2 to — Xi and — ® 2 . 
When we add the two loop diagrams, then terms involving odd total powers 
of N + p + q cancel, just as they should. 


18 



For other groups the diagrammatics is more complicated (see [l]) and 
in general the nonorientable diagram is essential for unitarity. We will not 
pursue this here - the group theoretic factors should not qualitatively alter 
our result. 

As is well known, the nonplanar orientable diagram we have included 
has singularities corresponding to propagation of closed string intermediate 
states (as may be seen heuristically by lifting the inner circle of Figure 5b). 
In fact it was this observation, and the requirement that the singularities be 
poles that led to the realisation of the importance of d = 26 in string theory. 
At tree level, this corresponds to a two point coupling between open and 
closed string states, necessary for unitarity, which means that open string 
states can in principle mix with closed string states to form the true mass 
eigenstates of the theory. However in our case it is clear that mixing cannot 
occur. This is because for the open string the states on the leading Regge 
trajectory have J = a M 2 + 1 and M 2 = ^r(AT — 1) whereas for the closed 
string the relation is J = |a*A f 2 + 2, and M 2 = A(AT — 1). It follows that 
for M 2 > 0 no closed string state has the correct mass and spin to mix 
with the leading Regge trajectory open string states. Thus we are not faced 
with singularities corresponding to intermediate closed string states and the 
consequent requirement that d = 26. 

5 Decay Rates 

Equation (26) with the added nonplanar diagram contains all the information 
pertaining to the decay rate of a leading trajectory state. However, it is not 
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in an immediately useable form for we have not given an analytic expression 
for the coefficients C^ q . To try to understand which terms dominate the sum 
in (26) we have written a computer program to explicitly compute the various 
coefficients for different external states, namely those with level number N 
up to 100 (Figure 6). We chose d=4 in order that the variation with N could 
clearly be seen. The results for N = 6 and N = 20 are shown in Figures 7 
and 8. The decay rate for a string at level N is dominated by decay into a 
photon (p, q = 1) and a massive string with q,p = N — 1. The rate for this 
process alone is also shown in Figure 6 as well as an analytic approximation 
to it we shall derive below. 

If we define, for given p, q 

& = ■£■£ />-<-■(»->) »(,) ( e 

r=0 twists n= — oo 

(30) 


where 


= F(JV + 1) 

' r(r + l)r(JV-r+l)r(f +r-i) 


(31) 


then the coefficient of x^x’ gives the rate of decay into a string at level p and 
one at level q. The crucial point is very simple. *P(r) is a rapidly decreasing 
function of r, and this effect dominates in d=4. 

Our strategy is to extract, for r = 0,1,2..., the coefficient of xf - 1 x 2 in 
(30). We denote this coefficient Q r . In this term, p = N — 1 and q = 1 or 
vica versa so from (25) a = 2 ^_ l y We find 


Q o 

Q i 


N 

(N- IV? 

- IN 

3 (jV — l)y/TT 
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jV(3A r - 4) 

24 (N - l) 2 v^ 

—N(N - 2)(19JV — 40) 
1260(JV - l) 3 ^ 


The sum of these terms appears to converge rapidly with increasing r 
The terms above add to 0.443/^ for large N. so that from (26) we have 

1 0.443<7 2 


IV- 


>N- 1,1 ~ 


(33) 


E ext 167T 

including an extra factor of two because of decay into x 2 x 2 and 
a factor of two from the optical theorem and a ‘symmetry’ factor of | (see 
Appendix B). This result is compared with the exact result as evaluated on 
a computer in Figure 6. 

This result certainly provides a lower bound on the total decay rate of 
this string state. The total decay rate is not dominated by this process for 
N < 100 . It does however appear to decrease more steeply as one goes 
to larger N - the magnitude of the slope of T versus N increases in the 


ratios 1 : 1.24 : 1.53 for the ranges 60 - 80,80 - 90,90 - 100, faster than 
a simple power law (the magnitude of the slope any power law decreases 
with increasing N). If we extrapolate the slope at N = 100 to higher N 
we would find the total decay rate reached (33) at N ~ 270. The total 
decay rate certainly must flatten out for N of this order. Unfortunately 
the exact computation becomes prohibitive at such large N . It does seem 
possible therefore that the total decay rate approaches (33) (possibly times a 
constant) for very large N . It seems very difficult to check this conjecture in 
our approach. In fact Green and Veneziano [5] argued that the total decay 
rate should go like l/E ex t - our result seems consistent with this. However 
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they applied this to meson resonances with N < 5 - we have shown a far 
better way to calculate the result for small N. 

In any case our results are rather surprising from the point of view of 
string interactions as splitting and joining - our result certainly conflicts 
with the idea of a splitting probability per unit length. This would lead to 
T oc y/N . Unless something very strange happens at large N we have the 
bound const/ \/N < T < const . 

If the lower bound is saturated for very large N this means that the decay 
is completely dominated by the ends of the string - if one determines from 
(33) a rate of energy loss 


P = TE. 


ext 


0.443$ 2 

167T 


(34) 


then this is completely independent of the mass (and thus the *size’[2] of the 
rotating string). Thus the lifetime of a string would be proportional to its 
length. 

It would also be interesting to calculate decay rates for strings in more 
general states, to see whether and how the picture changes. 

In Figures 7 and 8 we show plots of the amplitudes of the various terms 
in (30) for IV = 6 and N = 20. 

As we discussed in Section 1, we expect the decay rate to be very much 
suppressed in higher dimensions from phase space. We calculated the total 
decay rate in d = 8 and d = 26 for N = 20 and 40 to check this. The result, 
in the same units as Figure 6 but with appropriate powers of s/2a’ to m ak e 
up dimensions, were for d = 8: r 20 = 1.37. 10" 3 , r 40 = 1.31. 10" 3 and for 
d = 26: r 20 = 3.78.10 -lo ,r 4O = 7.79.10 -10 , much smaller as expected. These 
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cases are also different in that the N — 1, 1 term that we saw was largest in 
d = 4 is no longer largest here. These results give some indication that in 
high dimensions string cosmology may be rather complicated - large strings 
can have quite long lifetimes. Of course all this would also be affected by 
powers of the compactification radius in any ‘realistic’ model of open strings. 


6 Superstrings 


In this section we sketch the calculation for the decay rate of the analagous 
first excited level in the open superstring model. This is harder because 
the general vertex is not known, and the fermionic traces difficult. Thus we 
shall only deal with the simplest case. There are no surprises here - the first 
massive state decays in a Planck time. As in the bosonic case we use the 
optical theorem, and thus we compute the imaginary part of the two point 
one-loop amplitude. We perform the calculation of the amplitude using the 
vertex operator formalism in the light-cone gauge. As a first step we have 
to determine the vertex, in the light-cone gauge, for the emission of the first 
excited state of mass m = s/l (we set 2a = 1 throughout this section). To 
do this we write the most general combination of P* and operators [l] 

with conformal spin 1 compatible with the symmetric polarisation tensor 

V(C\ k l ) = gC j (\iP i P j + X 3 (P i R il k l + P j R il k l ) + X 3 R il k l R’ n k n )e ikX (35) 
with 


P* = dX'/dr 

R ij = (36) 
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Our notation as before is that of [1], is the traceless symmetric polari- 
sation tensor with with ( x ^k^ =■ = 0. The coefficients A 1? A 2 and A 3 can 
be determined by evaluating the matrix element of the coupling between two 
massless particles and one massive particle in two different ways. First we 
take the above vertex between two massless states 

G <h,n\V(C!,k[)\Km> & 

= sMCjCj + »A>[(*.<j)*>fl’« + (*>C0GCi’^ - (*.Ca)*5Ci - 

- yAsfcJGCi’G (37) 

(we used momentum conservation = 0) or by taking the matrix element 
of the vertex for the emission of a massless state between an excited and a 
massless state in a cyclic permutation of (37) 

C 2 C 3 < ki,n | g(P i - R'>ki)e' k > x a‘_ 1 | k u m > (? l 
= gCitfG + steMcM - gik&KitfH + g(GG)kid j ki (38) 
Comparing the results of (37) and (38) we find A x = 1,A 2 = A and A 3 = A. 
With the above vertex we can compute the two point one-loop amplitude 
defined as 

A k = Tr( AjVfcAaV.*) (39) 

where V k is the vertex 35 for the emission of a particle corresponding to 
the first excited state and A„ is the string propagator between neighboring 
vertices 

A a = \ / dx a x L a ~ x (40) 

z Jo 

and 

L = If 2 + £ x »< + (41) 

L n=l 
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Due to the fact that TrR n = 0 for n < 3 only the last term, proportional 
to A 3 contributes to the amplitude. Due to supersymmetry the partition 
function, f(w ) = n£=i(l _ ™ n ), cancels out. The trace over the oscillators 
can be easily carried out. 

For the amplitude we get 

Ak = [ Jf 10 P fix. f C 

16 J Jo Jo 

oo 

c = n 


.(1 — w 


n—1 


x 2 )(l - W n JX2 ) |2 


„=i (l-u; n ) 2 

using k 2 = —2 with w = x.x 2 and the kinematical factor 


(42) 

(43) 


K = C’( lTn k n k q k t k t & itl<imt 


(44) 


where 


t injti qmz _ Tr[R^Ri l R lq R^ z ) (45) 

as defined in [1]. Using the mass shell and gauge conditions we can evaluate 

K, 

K = 8C 2 = VS (46) 

using ( 2 — JE~~ = ^8 ^he momentum integration is carried out as in section 
2 and the polynomial C expanded as in section 3. The only term contributing 
an imaginary part is the first; it gives 1. We thus get 


i s’tf -3 !' r 


f dy f°° dx x- 4 e- ixA ~‘ 
Jo J-iK 


(2irp 16 * “• " (4?) 

with A = y 2 —y. The imaginary part is then, following the method presented 
in section 1 


1 irg 2 K f 1 

1 (27t) 5 16.6 


/ dy{y - y 2 ) 3 = i 
Jo 


9 2 Ktt 


(2ir) 5 16.6.140 


(48) 
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The same result is obtained from the diagram with one twist (which has an 
additional minus sign) , and the diagram with two twists. 

The first excited level decays into 2 massless ground state particles. The 
decay rate, in units where 2a — 1, is given by 


g 2 tt-v/2 

(2tt) 5 6720 


(49) 


The extension of the above calculation to higher excited states appears to be 
difficult in this formalism - being complicated in particular by the fermionic 
traces. 


7 Appendix A 

In general the function A(y ) is a parabola with zeros at y±. Shifting y to 
y = y — and defining a = we obtain 

In = f dy ( a 2 - y*) n 
J —a 

Changing variables to z = (j() 2 this is 

/„ = a 2n+1 f 1 dzz~*(l - z) n 
Jo 

which is just 

7 _2n+l TDf 1 _ , _an+l lr * r ( n + 1 ) 

U-a B(-,n + l) = c r(n + §) 
the required result. 


26 



Oft/QWAL 

0F POOR 


0 *or f S 

QUALITY 


8 Appendix B 

As we argued in Section (4) the two one loop diagrams must be added with 
equal weight. The overall weight of these diagrams can be found by explic- 
itly calculating the tree amplitudes for the lowest excited states and then 
comparing this answer to that derived via the optical theorem. We shall 
choose the simplest well-defined case - the decay of an N = 2 state into 
photons. The result we find is that if the tree level coupling constant is g 
(| for each diagram, see section 4) then each of the loop diagrams must be 
added in with a coupling (§) 2 i.e. one could take one loop diagram with a 
symmetry factor of | and the rule that C be conserved. If nothing else this 
appendix should convince the reader that calculating total decay rates by 
summing over polarisations and integrating over phase space is considerably 
more complicated than using the optical theorem! We shall also calculate in 
arbitrary d as a useful check of the formalism above. By calculating , C\ x 
and explicitly we find that 

r M> , = ^ET- PO* +!)(<*- 1) - !«(<* + !) + (<* + 6)1 f( 2 Ej ( 50) 

where an extra factor of has been added to correctly normalise the 

external state. Here is the polarisation tensor of the initial massive state. 

Now we check this using the tree amplitude method. The polarisation 
vectors of the photons will be called e a , e 7 . The matrix element for the 
transition is given by 

T kX = g < k 3 | a“ : P*PV* J * : al | k 2 ~ k i ~ **) (51) 

It is important to remember all the constraints on the polarisation tensor and 
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vectors, namely that they are transverse and that the tensor has no trace. 
Calculating all the various terms in (51) we arrive at 

T kX = g(v ai, ^k x + 2t7 q V 7 ~ 2k 3 aT)^k x + 2k^ v aK k x - k^k x k^k^) 

6a4 c ^(*3 + k x - k 2 ) (52) 

and so the decay rate is simply 

r = 7 fE/ - h - fc 3 ) | T~ x | ! (53) 

2 C 1 ^ 3 

where we sum over the final polarisation vectors and include a factor of | for 
the two identical final particles. The various identities we will need are: 


kU 




E 


= 0 


^ cl 6l 


k a k & 


(54) 


where the arrow over any quantity indicates that only the spatial components 
are nonzero. We now choose the frame k 2 = (M, 0, ...) which implies that 
k 2 = 0 so 


*?E<-4 = 0 

€ 

Hence we can rewrite (53) 
r - « 

(2j? a V 7 + v^Kk^TT-o* + v*kikl)(6^ - 

*i 

- 2 JvcAC/tA 


(55) 


(56) 
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but E 2 = V2, thus after a little algebra we find 


= _L_ f 

l6x/2 J 


16v/2 J El 


f(E 2 - 2f?j) [4 ij*V a - 


* AA 


AA 


fc l 


-4 k^rf^ 

kl 

-Ak^klkl 

kl 


kl 

+ {d — 2)k*k* kl £ k % 


we now use the identity 


(57) 


2n indices _ 

I*"* = (27r) d - 2 r(^^) 

where “spare” indices from the LHS are used up 
of metric tensors e.g 


J k 2n+d ~ 2 f(k 2 )dk (58) 

in all distinct permutations 


J d^kl klklklkif{k 2 ) 


(2*)*-*r(^) 

J k d+2 f(kl ) (ifV* + ^V*) (59) 


all other 7 ]T] giving zero due to the tracelessness of Doing this for all 

terms we recover exactly (50). 
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Figure Captions 

Figure 1 

The optical theorem in lowest order for a coupling. 

Figure 2 

Contours for integration in determining the imaginary part of the one loop 
amplitude. 

Figure 3 

The function A(y) 

Figure 4 

Constructing the planar self energy loop by sewing the ends of a tree 
diagram. 

Figure 5 

(a) The nonorientable nonplanar self energy diagram (b) The orientable 
nonplanar self energy diagram Diagram (a) is excluded for U(l) open 
strings since cutting it reveals a string with two like charges on the ends. 

Figure 6 

The total decay rate in d = 4 as a function of level number N up to 
N = 100 (crosses on solid line). It is given in units of (g 2 / I6ic)y/2a' . Also 
shown is the decay rate into a single photon and a state at level N — 1, the 
largest single process, in the same units (crosses on dashed line), and our 
analytic approximation to it, equation (33) (plain dashed line), which 
becomes very good for large N. 

Figure 7 

The decay of an N = 6 string, p and q correspond to the level numbers of 
the decay products - the axes are marked off in integers starting at 
p = q = 1, corresponding to massless particles. p + q = odd decays are 
forbidden by C conservation. The vertical scale is in units of the height of 
the largest peaks, which correspond to decay into a photon plus a string at 
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the next level down. 

Figure 8 

The decay of an N = 20 string. The axes are as in Figure 7. 
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